Where Are the Formerly Y-linked Genes in the Ryukyu Spiny Rat that has Lost its Y Chromosome?

It has been predicted that the highly degenerate mammalian Y chromosome will be lost eventually. Indeed, Y was lost in the Ryukyu spiny rat Tokudaia osimensis, but the fate of the formerly Y-linked genes is not completely known. We looked for all 12 ancestrally Y-linked genes in a draft T. osimensis genome sequence. Zfy1, Zfy2, Kdm5d, Eif2s3y, Usp9y, Uty, and Ddx3y are putatively functional and are now located on the X chromosome, whereas Rbmy, Uba1y, Ssty1, Ssty2, and Sry are missing or pseudogenized. Tissue expressions of the mouse orthologs of the retained genes are significantly broader/higher than those of the lost genes, suggesting that the destinies of the formerly Y-linked genes are related to their original expressions. Interestingly, patterns of gene retention/loss are significantly more similar than by chance across four rodent lineages where Y has been independently lost, indicating a level of certainty in the fate of Y-linked genes even when the chromosome is gone.

[Comparison on the diversity of antibiotic resistance genes of three rodent species]. / 3ç§å•®é½¿åŠ¨ç‰©æŠ—ç”Ÿç´ æŠ—æ€§åŸºå› å¤šæ ·æ€§æ¯”è¾ƒ.

Antibiotic resistance genes (ARGs) have attracted widespread attention as a new global pollutant, mainly due to the abuse of antibiotics. To investigate the diversity of ARGs in three rodent species, we used metagenomic sequencing analysis to analyze the diversity of antibiotic resistance genes of 17 individuals of Apodemus peninsulae and 17 individuals of Myodes rufocanus collected from Mudanfeng, and nine individuals of Apodemus agrarius collected from Sandaoguan. A total of 19 types and 248 subclasses of ARGs were detected in the three rodent species. Seven ARGs showed significant difference and five ARGs showed extremely significant difference between M. rufocanus and A. agrarius. Seven ARGs showed significant difference and four ARGs showed extremely significant difference between A. peninsulae and A. agrarius. Four ARGs showed significant difference and five ARGs showed extremely significant difference between M. rufocanus and A. peninsulae. ARGs showing high abundance in three rodents were macrolides, lincoamides, tetracyclines, and ß-lactams. ARGs were widely distributed in the three rodent species. The significant differences in ARGs among different species might be due to the different distribution areas and their diet differentiation. The study could provide a basis for further studies of ARGs in mice and improve the understanding of the harm of ARGs transmission.

Speciation across biomes: Rapid diversification with reproductive isolation in the Australian delicate mice.

Phylogeographic studies of continental clades, especially when combined with palaeoclimate modelling, provide powerful insight into how environment drives speciation across climatic contexts. Australia, a continent characterized by disparate modern biomes and dynamic climate change, provides diverse opportunity to reconstruct the impact of past and present environments on diversification. Here, we use genomic-scale data (1310 exons and whole mitogenomes from 111 samples) to investigate Pleistocene diversification, cryptic diversity, and secondary contact in the Australian delicate mice (Hydromyini: Pseudomys), a recent radiation spanning almost all Australian environments. Across northern Australia, we find no evidence for introgression between cryptic lineages within Pseudomys delicatulus sensu lato, with palaeoclimate models supporting contraction and expansion of suitable habitat since the last glacial maximum. Despite multiple contact zones, we also find little evidence of introgression at a continental scale, with the exception of a potential hybrid zone in the mesic biome. In the arid zone, combined insights from genetic data and palaeomodels support a recent expansion in the arid specialist P. hermannsburgensis and contraction in the semi-arid P. bolami. In the face of repeated secondary contact, differences in sperm morphology and chromosomal rearrangements are potential mechanisms that maintain species boundaries in these recently diverged species. Additionally, we describe the western delicate mouse as a new species and recommend taxonomic reinstatement of the eastern delicate mouse. Overall, we show that speciation in an evolutionarily young and widespread clade has been driven by environmental change, and potentially maintained by divergence in reproductive morphology and chromosome rearrangements.

A recent gibbon ape leukemia virus germline integration in a rodent from New Guinea.

Germline colonization by retroviruses results in the formation of endogenous retroviruses (ERVs). Most colonization's occurred millions of years ago. However, in the Australo-Papuan region (Australia and New Guinea), several recent germline colonization events have been discovered. The Wallace Line separates much of Southeast Asia from the Australo-Papuan region restricting faunal and pathogen dispersion. West of the Wallace Line, gibbon ape leukemia viruses (GALVs) have been isolated from captive gibbons. Two microbat species from China appear to have been infected naturally. East of Wallace's Line, the woolly monkey virus (a GALV) and the closely related koala retrovirus (KoRV) have been detected in eutherians and marsupials in the Australo-Papuan region, often vertically transmitted. The detected vertically transmitted GALV-like viruses in Australo-Papuan fauna compared to sporadic horizontal transmission in Southeast Asia and China suggest the GALV-KoRV clade originates in the former region and further models of early-stage genome colonization may be found. We screened 278 samples, seven bat and one rodent family endemic to the Australo-Papuan region and bat and rodent species found on both sides of the Wallace Line. We identified two rodents (Melomys) from Australia and Papua New Guinea and no bat species harboring GALV-like retroviruses. Melomys leucogaster from New Guinea harbored a genomically complete replication-competent retrovirus with a shared integration site among individuals. The integration was only present in some individuals of the species indicating this retrovirus is at the earliest stages of germline colonization of the Melomys genome, providing a new small wild mammal model of early-stage genome colonization.

Insulin and circadian rhythm genes of the Nile rat (Arvicanthis niloticus) are conserved and orthologous to those in the rat, mouse and human.

The African grass or Nile rat (NR) (Arvicanthis niloticus) is a herbivorous diurnal rodent which is used as a biological model for research on type 2 diabetes mellitus (T2DM) and the circadian rhythm. Similar to humans, male NRs develop T2DM with high-carbohydrate diets. The NR thus provides a unique opportunity to identify the nutritional and underlying genetic factors that characterise human T2DM, as well as the effects of potential anti-diabetic phytochemicals such as Water-Soluble Palm Fruit Extract. Whole genome sequencing (WGS) could help identify possible genetic causes why NRs spontaneously develop T2DM in captivity. In this study, we performed WGS on a hepatic deoxyribonucleic acid (DNA) sample isolated from a male NR using PacBio high-fidelity long-read sequencing. The WGS data obtained were then de novo assembled and annotated using PacBio HiFi isoform sequencing (Iso-Seq) data as well as previous Illumina RNA sequencing (RNA-Seq) data. Genes related to insulin and circadian rhythm pathways were present in the NR genome, similar to orthologues in the rat, mouse and human genomes. T2DM development in the NR is thus most likely not attributable to structural differences in these genes when compared to other biological models. Further studies are warranted to gain additional insights on the genetic-environmental factors which underlie the genetic permissiveness of NRs to develop T2DM.

Chromosomal-level assembly of Tokudaia osimensis, Tokudaia tokunoshimensis, and Tokudaia muenninki genomes.

Herein, we present the first high-quality long-read-based chromosome-level genome assemblies and gene annotations of the genomes of three endangered Tokudaia species: Tokudaia osimensis, Tokudaia tokunoshimensis, and Tokudaia muenninki. These species, which are endemic to different islands of the Ryukyu Islands, Japan, exhibited unique karyotypes and sex chromosomal characteristics. The genome assemblies generated using PacBio, Illumina, and Hi-C sequence data consisted of 13 (corresponded to 12 autosomes and one X chromosome), 23 (corresponded to 22 autosomes and one X chromosome), and 23 (corresponded to 21 autosomes and the neo- and ancestral X regions) chromosome-level scaffolds that contained 2,445, 2,477, and 2,661 Mbp of sequence data, respectively. Annotations of protein-coding genes were performed using RNA-Seq-based, homology-based, and Ab initio methods. BUSCO completeness values for every species exceeded 96% for genomes and 98% for genes. These data can be an important resource for contributing to our understanding of species genomes resulting from allopatric speciation and provide insights into mammalian sex-determination mechanisms and sex chromosome evolution.

Immunogenetics, sylvatic plague and its vectors: insights from the pathogen reservoir Mastomys natalensis in Tanzania.

Yersinia pestis is a historically important vector-borne pathogen causing plague in humans and other mammals. Contemporary zoonotic infections with Y. pestis still occur in sub-Saharan Africa, including Tanzania and Madagascar, but receive relatively little attention. Thus, the role of wildlife reservoirs in maintaining sylvatic plague and spillover risks to humans is largely unknown. The multimammate rodent Mastomys natalensis is the most abundant and widespread rodent in peri-domestic areas in Tanzania, where it plays a major role as a Y. pestis reservoir in endemic foci. Yet, how M. natalensis' immunogenetics contributes to the maintenance of plague has not been investigated to date. Here, we surveyed wild M. natalensis for Y. pestis vectors, i.e., fleas, and tested for the presence of antibodies against Y. pestis using enzyme-linked immunosorbent assays (ELISA) in areas known to be endemic or without previous records of Y. pestis in Tanzania. We characterized the allelic and functional (i.e., supertype) diversity of the major histocompatibility complex (MHC class II) of M. natalensis and investigated links to Y. pestis vectors and infections. We detected antibodies against Y. pestis in rodents inhabiting both endemic areas and areas considered non-endemic. Of the 111 nucleotide MHC alleles, only DRB*016 was associated with an increased infestation with the flea Xenopsylla. Surprisingly, we found no link between MHC alleles or supertypes and antibodies of Y. pestis. Our findings hint, however, at local adaptations towards Y. pestis vectors, an observation that more exhaustive sampling could unwind in the future.

Comparative genomics of widespread and narrow-range white-bellied rats in the Niviventer niviventer species complex sheds light on invasive rodent success.

Widespread species that inhabit diverse environments possess large population sizes and exhibit a high capacity for environmental adaptation, thus enabling range expansion. In contrast, narrow-range species are confined to restricted geographical areas and are ecologically adapted to narrow environmental conditions, thus limiting their ability to expand into novel environments. However, the genomic mechanisms underlying the differentiation between closely related species with varying distribution ranges remain poorly understood. The Niviventer niviventer species complex (NNSC), consisting of highly abundant wild rats in Southeast Asia and China, offers an excellent opportunity to investigate these questions due to the presence of both widespread and narrow-range species that are phylogenetically closely related. In the present study, we combined ecological niche modeling with phylogenetic analysis, which suggested that sister species cannot be both widespread and dominant within the same geographical region. Moreover, by assessing heterozygosity, linkage disequilibrium decay, and Tajima's D analysis, we found that widespread species exhibited higher genetic diversity than narrow-range species. In addition, by exploring the "genomic islands of speciation", we identified 13 genes in highly divergent regions that were shared by the two widespread species, distinguishing them from their narrow-range counterparts. Functional annotation analysis indicated that these genes are involved in nervous system development and regulation. The adaptive evolution of these genes likely played an important role in the speciation of these widespread species.

Spatial transcriptomics reveals asymmetric cellular responses to injury in the regenerating spiny mouse (Acomys) ear.

In contrast to other mammals, the spiny mouse (Acomys) regenerates skin and ear tissue, which includes hair follicles, glands, and cartilage, in a scar-free manner. Ear punch regeneration is asymmetric with only the proximal wound side participating in regeneration. Here, we show that cues originating from the proximal side are required for normal regeneration and use spatially resolved transcriptomics (tomo-seq) to understand the molecular and cellular events underlying this process. Analyzing gene expression across the ear and comparing expression modules between proximal and distal wound sides, we identify asymmetric gene expression patterns and pinpoint regenerative processes in space and time. Moreover, using a comparative approach with nonregenerative rodents (Mus, Meriones), we strengthen a hypothesis in which particularities in the injury-induced immune response may be one of the crucial determinants for why spiny mice regenerate whereas their relatives do not. Our data are available in SpinyMine, an easy-to-use and expandable web-based tool for exploring Acomys regeneration-associated gene expression.

Bartonella raoultii sp. nov., isolated from infected rodents (Mastomys erythroleucus) in Senegal.

Bartonella species are involved in various human diseases, causing a range of clinical manifestations; animals are considered as the main reservoirs, transmitting diverse species of Bartonella through direct contact and haematophagous insects. Here, we characterize a new species, Bartonella raoultii sp. nov., within the genus Bartonella, using a taxonogenomic polyphasic approach. Strain 094T (= CSUR B1097T=DSM 28004T), isolated from the blood of an infected rodent (Mastomys erythroleucus) in Senegal, is an aerobic and rod-shaped bacterium. The annotated non-contiguous genome sequence is 1 952322 bp long and contains 37.2 mol% G+C content, 1686 protein-coding genes and 50 RNA genes, including seven rRNA genes.

Genome Report: chromosome-scale genome assembly of the African spiny mouse (Acomys cahirinus).

There is increasing interest in the African spiny mouse (Acomys cahirinus) as a model organism because of its ability for regeneration of tissue after injury in skin, muscle, and internal organs such as the kidneys. A high-quality reference genome is needed to better understand these regenerative properties at the molecular level. Here, we present an improved reference genome for A. cahirinus generated from long Nanopore sequencing reads. We confirm the quality of our annotations using RNA sequencing data from 4 different tissues. Our genome is of higher contiguity and quality than previously reported genomes from this species and will facilitate ongoing efforts to better understand the regenerative properties of this organism.

Diversification and evolutionary history of the African laminated-toothed rats (Rodentia, Otomyini).

The African continent was subjected to periodic climatic shifts during the Pliocene and Pleistocene. These habitat changes greatly affected the evolutionary processes and tempo of diversification in numerous, widely distributed mammals. The Otomyini (Family Muridae) comprises three African rodent genera, Parotomys, Otomys and Myotomys, characterized by unique laminated-shaped molars. Species within this tribe generally prefer open-habitat and show low dispersal capabilities, with previous studies suggesting that their diversification was closely associated with climatic oscillations over the last four million years. Our phylogenetic reconstructions, based on three mitochondrial (mtDNA) genes (Cytb, COI and 12S) and four nuclear introns (EF, SPTBN, MGF and THY), identified eight major genetic clades that are distributed across southern, eastern and western Africa. Our data permit the re-examination of the taxonomic status of the three genera as well as the previously proposed mesic-arid dichotomy of the 10 South African species. Moreover, multiple mtDNA species delimitation methods incorporating 168 specimens estimated the number of Otomyini species to be substantially higher than the âˆ¼ 30 recognized, suggesting that the current taxonomy will necessitate an integrative approach to delimit extant species diversity within the Otomyini. The data suggests that the origin of the tribe can be dated back to âˆ¼ 5.7 million years ago (Ma) in southern Africa. The distribution and phylogenetic associations among the eight major otomyine evolutionary lineages can best be explained by several waves of northward colonization from southern Africa, complemented by independent reversed dispersals from eastern back to southern Africa at different time periods. There is strong support for the hypothesis that the radiation, dispersion, and diversification of the otomyine rodents is closely linked to recent Plio-Pleistocene climatic oscillations.

The effects of antibiotics and illness on gut microbial composition in the fawn-footed mosaic-tailed rat (Melomys cervinipes).

The gut microbiota are critical for maintaining the health and physiological function of individuals. However, illness and treatment with antibiotics can disrupt bacterial community composition, the consequences of which are largely unknown in wild animals. In this study, we described and quantified the changes in bacterial community composition in response to illness and treatment with antibiotics in a native Australian rodent, the fawn-footed mosaic-tailed rat (Melomys cervinipes). We collected faecal samples during an undiagnosed illness outbreak in a captive colony of animals, and again at least one year later, and quantified the microbiome at each time point using 16s ribosomal rRNA gene sequencing. Gut bacterial composition was quantified at different taxonomic levels, up to family. Gut bacterial composition changed between time periods, indicating that illness, treatment with antibiotics, or a combination affects bacterial communities. While some bacterial groups increased in abundance, others decreased, suggesting differential effects and possible co-adapted and synergistic interactions. Our findings provide a greater understanding of the dynamic nature of the gut microbiome of a native Australian rodent species and provides insights into the management and ethical well-being of animals kept under captive conditions.

Circadian functioning of Locus Cœruleus of the nocturnal rat and diurnal rodent Arvicanthis.

The noradrenergic Locus Cœruleus is one of the major arousal structures involved in inducing wakefulness. While brain noradrenaline (NA) amounts display 24-h variations, the origin of NA rhythm is currently unknown. In this study, we tested the hypothesis that NA rhythm could result from its rhythmic synthesis. Therefore, we investigated the 24-h expression profile of NA rate-limiting enzyme, tyrosine hydroxylase (th), in the Locus Cœruleus (LC) of the nocturnal rat and the diurnal rodent Arvicanthis, under 12 h:12 h light/dark (LD) and constant darkness (DD) conditions. In both species, th mRNA levels vary significantly over 24-h. In nocturnal rats, th mRNA profiles show a unimodal rhythm, with peak values in late day in LD, and in the middle of the subjective day in DD. In contrast, th mRNA rhythm in Arvicanthis is characterized by a bimodal profile, with higher levels at the beginning of the day and of the night in LD, and in the middle of the subjective day and night in DD. The rhythmic pattern of th expression may be dependent on a LC clock machinery. Therefore, we investigated the expression of three clock genes, namely bmal1, per1, and per2, and found that their mRNAs display significant variations between day and nighttime points in both species, but in opposite directions. These data show that NA rhythm may be related to circadian expression of th gene in both species, but differs between nocturnal and diurnal rodents. Furthermore, the phase opposition of clock gene expression in the rat compared to Arvicanthis suggests that the clock machinery might be one of the mechanisms involved in th rhythmic expression.

Distribution of vasopressin 1a and oxytocin receptor protein and mRNA in the basal forebrain and midbrain of the spiny mouse (Acomys cahirinus).

The nonapeptide system modulates numerous social behaviors through oxytocin and vasopressin activation of the oxytocin receptor (OXTR) and vasopressin receptor (AVPR1A) in the brain. OXTRs and AVPR1As are widely distributed throughout the brain and binding densities exhibit substantial variation within and across species. Although OXTR and AVPR1A binding distributions have been mapped for several rodents, this system has yet to be characterized in the spiny mouse (Acomys cahirinus). Here we conducted receptor autoradiography and in situ hybridization to map distributions of OXTR and AVPR1A binding and Oxtr and Avpr1a mRNA expression throughout the basal forebrain and midbrain of male and female spiny mice. We found that nonapeptide receptor mRNA is diffuse throughout the forebrain and midbrain and does not always align with OXTR and AVPR1A binding. Analyses of sex differences in brain regions involved in social behavior and reward revealed that males exhibit higher OXTR binding densities in the lateral septum, bed nucleus of the stria terminalis, and anterior hypothalamus. However, no association with gonadal sex was observed for AVPR1A binding. Hierarchical clustering analysis further revealed that co-expression patterns of OXTR and AVPR1A binding across brain regions involved in social behavior and reward differ between males and females. These findings provide mapping distributions and sex differences in nonapeptide receptors in spiny mice. Spiny mice are an excellent organism for studying grouping behaviors such as cooperation and prosociality, and the nonapeptide receptor mapping here can inform the study of nonapeptide-mediated behavior in a highly social, large group-living rodent.

Phylogeography of the striped field mouse (Apodemus agrarius Pallas, 1771) in light of new data from central part of Northern Eurasia.

A phylogeographic analysis of A. agrarius based on the complete mtDNA cytochrome b and control region sequences has been performed using data obtained for the first time for the species from large regions of the central part of Northern Eurasia (23 localities of Altai, Western Siberia, and the Urals). The obtained results have demonstrated a complex intraspecific differentiation of A. agrarius, which has manifested not only in the isolation of the isles populations in Southeast Asia (Jeju and Taiwan), but also in the genetic heterogeneity of mainland populations, which has reflected the history of the modern intraspecific genetic diversity formation against the background of changing physiographic conditions of Eurasia in the Quaternary. The divergence of genetic lineages has taken place apparently simultaneously (in mid-Pleistocene) on the territory of the Eastern part of the modern disjunctive range, where all the identified lineages are present today. The demographic history and possible evolutionary scenarios for A. agrarius in the Western part of the range have been considered. TMRC reconstructions have shown that the lifetime of the common ancestor of the lineage that expanded in the Western Palearctic is about 17.7 [95% HPD 13.2-22.5] kyr. This suggests that the transcontinental expansion of A. agrarius is a relatively recent event that has occurred after the LGM.

Chromosomal Aberrations in Large Japanese Field Mice (Apodemus speciosus) Captured in Various Periods after Fukushima Daiichi Nuclear Power Plant Accident.

After the Fukushima Daiichi Nuclear Power Plant accident, we studied the chromosomal aberrations (dicentrics and translocations) in the splenic lymphocytes of wild mice inhabiting Fukushima prefecture. Here, we report the frequencies of chromosomal aberrations in large Japanese field mice (Apodemus speciosus) captured from 2012 to 2016 in a heavily contaminated area. The chromosomal aberrations were detected using newly developed 4-color FISH (fluorescence in situ hybridization) with A. speciosus chromosome 1-, 3-, 4- and 5-specific painting probes. The frequencies of chromosomal aberrations in mice captured in July 2012 and October 2014 were significantly higher than that in the mice inhabiting the non-contaminated control area; however, the frequency of chromosomal aberrations in mice captured in January 2016 was not. The frequency of chromosomal aberrations in individual mice tended to increase with certain dose rates and accumulated doses. Regression tree analyses suggested increasing chromosomal aberration rate in mice exposed to chronic radiation at dose rates of more than 1.1 mGy day-1 and at accumulated doses of more than 200 mGy. It is concluded that ambient dose rates in the most severely contaminated area of Fukushima prefecture and radiation doses to wild mice inhabiting this area decrease with time; consequently, chromosomal aberrations induced by radiation have not been detected 5 years after the accident.

Intrinsic differences in rod and cone membrane composition: implications for cone degeneration.

PURPOSE: In many retinal pathological conditions, rod and cone degeneration differs. For example, the early-onset maculopathy Stargardts disease type 1 (STGD1) is typified by loss of cones while rods are often less affected. We wanted to examine whether there exist intrinsic membrane differences between rods and cones that might explain such features. METHODS: Abca4 mRNA and protein levels were quantified in rod- and cone-enriched samples from wild-type and Nrl-/- mice retinas; rod- and cone-enriched outer segments (ROS and COS respectively) were prepared from pig retinas, and total lipids were analyzed by flame ionization, chromatography, and tandem mass spectrometry. Immunohistochemical staining of cone-rich rodent Arvicanthis ansorgei retinas was conducted, and ultra-high performance liquid chromatography of lipid species in porcine ROS and COS was performed. RESULTS: Abca4 mRNA and Abca4 protein content was significantly higher (50-300%) in cone compared to rod-enriched samples. ROS and COS displayed dramatic differences in several lipids, including very long chain poly-unsaturated fatty acids (VLC-PUFAs), especially docosahexaenoic acid (DHA, 22:6n-3): ROS 20.6% DHA, COS 3.3% (p < 0.001). VLC-PUFAs (> 50 total carbons) were virtually absent from COS. COS were impoverished (> 6× less) in phosphatidylethanolamine compared to ROS. ELOVL4 ("ELOngation of Very Long chain fatty acids 4") antibody labelled Arvicanthis cones only very weakly compared to rods. Finally, there were large amounts (905 a.u.) of the bisretinoid A2PE in ROS, whereas it was much lower (121 a.u., ~ 7.5-fold less) in COS fractions. In contrast, COS contained fivefold higher amounts of all-trans-retinal dimer (115 a.u. compared to 22 a.u. in rods). CONCLUSIONS: Compared to rods, cones expressed higher levels of Abca4 mRNA and Abca4 protein, were highly impoverished in PUFA (especially DHA) and phosphatidylethanolamine, and contained significant amounts of all-trans-retinal dimer. Based on these and other data, we propose that in contrast to rods, cones are preferentially vulnerable to stress and may die through direct cellular toxicity in pathologies such as STGD1.

[Analysis of intestinal microbial diversity in Leopoldamys edwardsi based on illumina sequencing technique].

To explore the composition and diversity of the intestinal microflora of Leopoldamys edwardsi in Hainan Island. In November 2019, DNA was extracted from fecal samples of 25 adult Leopoldamys edwardsi (14 males and 11 females) in Hainan Island at the Joint Laboratory of tropical infectious diseases of Hainan Medical College and Hong Kong University. Based on the IonS5TMXL sequencing platform, single-end sequencing (Single-End) was used to construct a small fragment library for single-end sequencing. Based on Reads shear filtration and OTUs clustering. The species annotation and abundance analysis of OTUs were carried out by using mothur method and SSUrRNA database, and further conducted α diversity and ß diversity analysis. A total of 1481842 high quality sequences, belonging to 14 Phyla, 85 families and 186 Genera, were obtained from 25 intestinal excrement samples of Leopoldamys edwardsi. At the level of phyla classification, the main core biota of the Leopoldamys edwardsi contained Firmicutes (46.04%),Bacteroidetes (25.34%), Proteobacteria (17.09%), Tenericutes (7.38%) and Actinobacteria (1.67%), these five phyla account for 97.52% of all phyla. The ratio of Helicobacter which occupied the largest proportion at the genus level was 12.44%, followed by Lactobacillus (11.39%), Clostridium (6.19%),Mycoplasma (4.23%) and Flavonifractor (3.52%). High throughput sequencing analysis showed that the intestinal flora of Leopoldamys edwardsi in Hainan Island was complex and diverse, which had the significance of further research.

Sympatric speciation of the spiny mouse from Evolution Canyon in Israel substantiated genomically and methylomically.

SignificanceWhether sympatric speciation (SS) is rare or common is still debated. Two populations of the spiny mouse, Acomys cahirinus, from Evolution Canyon I (EC I) in Israel have been depicted earlier as speciating sympatrically by molecular markers and transcriptome. Here, we investigated SS both genomically and methylomically, demonstrating that the opposite populations of spiny mice are sister taxa and split from the common ancestor around 20,000 years ago without an allopatric history. Mate choice, olfactory receptors, and speciation genes contributed to prezygotic/postzygotic reproductive isolation. The two populations showed different methylation patterns, facilitating adaptation to their local environment. They cope with abiotic and biotic stresses, due to high solar interslope radiation differences. We conclude that our new genomic and methylomic data substantiated SS.